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I. INTRODUCTION 
Hydrodynamicequationsinclude continuityequationand 
Euler’s momentum equation formulated using the Taylor 
series O(51) (Dean (1991). Meanwhile, KFSBC is a total 
velocity equation of the movement of the water surface in 
the direction of vertical axis that can be formulated using 
the Taylor series. 


Analytical solutions to Laplace’s equation using the 
separation of variables produce a sinusoidal wave equation 
(Dean (1991). Thus, the formulation of the equations for 
the water wave mechanics should be based on the nature of 
the sinusoidal function. Time step and spatial interval in 
the Taylor series for sinusoidal equations are correlated 
with phase speed (Courant (1928)), not with particle 
velocity. Hence, it is necessary to formulate a Taylor series 
in which time step and spatial intervals correlate with the 
water particle velocity. Thus, it can be used in the 
formulation of basic equations of hydrodynamics that are 
the basic equations of water wave mechanics. 


The first step of this research was formulating the Taylor 
series for sinusoidal functions where the time step and 
spatial interval can be correlated with the water particle 
velocity. At this stage, the weighted Taylor series was 
produced, that is, the Taylor series in which, there is a 
weighted coefficient on the time step and spatial interval. 


www.ijaers.com 


Abstract— In this study, the Taylor series is formulated with a weighted 
coefficient to time step and spatial interval. With the weighted Taylor 
series, the weighted total acceleration is formulated on Euler’s momentum 
equation and the Kinematic Free Surface Boundary Condition (KF SBC). 


The final part is the development of a time series water wave model using 
the weighted momentum Euler equation and the weighted KFSBC. 


Next, with the weighted Taylor series, the basic 
equations of hydrodynamics were formulated. They are 
namely the continuity equation, the Euler’s momentum 
equation,and KFSBC containing a weighted coefficient. 


With the basic equations of hydrodynamicscontaining the 
weighted coefficient, the time series water wave model 
was formulated. 


Il. THE FORMULATION OF THE WEIGHTED 
TAYLOR SERIES 


: ae: : ôx. : : 
This chapter examining the meaning of yin a sinusoidal 


f : 6z. : . 
function and the meaning of sin the hyperbolic functions 


considering the solution of Laplace’s equation which is the 
multiplication of a sinusoidal function with a hyperbolic 
function (Dean (1991). This chapter is a rewrite of 
Hutahaean (2021), considering that this section is the basis 
of the theory developed and at the same time is a 
correction of typos in Hutahaean (2021). 


2.1. An Overview of the Solution of Laplace's equation 
Solution of Laplace’s equation (Dean (1991) is, 
@(x, z,t) = Gcoshk(h + z)coskxsinot.....(1) 


Particle velocity in the direction of horizontal axis—xis, 


“= z2 = Gkcoshk(h + z)sinkxsinot......(2) 
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The velocity in the direction of vertical axis-zis 


w=- = = —Gksinhk(h + z)coskxsinot...(3) 


k : wave number = = (m!) 

L : wavelength (m) 

o : angular frequency = = (sec’! ) 
T : wave period (sec.) 

h : water depth (m) 


i ; 8x. F 
From Laplace’s equation, = in the Taylor series for a 


sinusoidal water wave equation is not the water particle 
velocity, it should be the wave celerity or wave phase 


.). OZ. ; f 
speed. Meanwhile Sis also a function of wave celerity that 


is described in the following section. 


2.2. A function of a single variable 


The first step was examining the characteristics of ôt, xin 
the sinusoidal function and ôzin the hyperbolic function in 
the Taylor series, in a function of a single variable.The 
formula of the Taylor series for a function with one 
variable is: 


PEE EREE E a 4 AR eee 





dx 2! dx? 3! dx? 
éx* atf ôx” af 
ggat eee t r gn e 
a. fŒ) =cosot 


The first single-variable of sinusoidal function 
examinedwas f (t) = cosi/olot. In his function, the value of 
dtwas examined, in which the Taylor series can be used 
with only one derivative. This study was carried out using 
the Taylor series third order, 

St? d?f | dt d?f 

2! dt? 3! ats" 





f(t + dt) =f) + st 4+ (5) 


The second and third differential terms can be ignored if 

the sum of the two terms is much smaller than the first 

term: 

strap, st3af 

2! dt? 3! dt3 
ace 








< ¢.....(6) 








The fourth term, fifth term, and so on can be used. 
However, considering that ôt is a very small number, the 
fourth and higher differential term is a very small number 
that can be ignored. Equation (6) is hereinafter referred to 
as the optimization equation. In (6), the variable to be 
calculated is dt. While € is a very small number which will 
determine the level of accuracy. dt in the denominator 
with the numerator cancel each other out, 
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Std2f dst*d3f 
2!dt? 3! dt? 
df 








<E 


dt 








The derivatives of the function are 


af 


: a? a3 
—=-oSsinot ; ai f 
dt 2 


= —o’cosot dan — = o?sinot. 
t dt? 
The substitution of the derivative of the function in (7), 


2 
2 (—0o°cosot) + 2 (o°sinot) 
2 6 <e 





—osinot 


This equation is valid for any value of ot as long as it is 
not equal to zero. It is easier to use the value of otwhere 
sinot = cosot.This is called the characteristic point. The 
final equation is: 


ôt jot 
o> =g — 


< 
2 ei 








For very small dt, the term in the absolute value sign will 
be positive. Thus, the absolute sign can be omitted, 


By using an equal sign, 
2 
-2 624+ 26t—-2=0 2 on (8) 
6 2 
Equation (8) is for calculating ôt where the Taylor series 
can be used only with the first differential. 
b. f(x) = cos kx 
Next, dxwas calculated in the functionf (x) = cos kx. In 
the same way,the obtained formula is, 
2 
— © 8x? + £6x — £ = 0....9) 
6 2 
Equation (9) is for calculating 6x where the Taylor series 
can be used only with the first differential. 
c. f(z) = coshk(h + z) 
The function of the next variable isf (z) = cosh k(h + 2). 
In the same way,the obtained formula is, 
2 
F622 +É58z-— e = 0.....(10) 
6 2 
With (10),6z, can be calculated, where the Taylor series 
can be used only with the first differential. 


In Table (1), it is presented the calculation result ofdt, dx, 
and 6z, with (8), (9), and (10), in which wave number 
kcalculated using the dispersion equation of the linear 
wave theory, at waterdepth of h =10 m.The dispersion 
equation of the linear wave theory (Dean (1991) is, 


o*=gktanhkh _...... (11) 


g ‘gravitational force 
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Table.1: The calculation results ofSt, 5x, and 6z 



































T ôt ôx Oz 
(sec.) (sec.) (m) (m) 
6 0,00191 | 0,01542 | 0,0154 
7 0,00223 | 0,01905 | 0,01903 
8 0,00255 | 0,02258 | 0,02255 
9 0,00287 | 0,02603 0,026 
10 0,00319 | 0,02942 | 0,02938 
11 0,0035 | 0,03277 | 0,03273 
12 0,00382 | 0,03609 | 0,03604 
13 0,00414 | 0,03938 | 0,03933 
14 0,00446 | 0,04265 | 0,04259 
15 0,00478 | 0,04591 | 0,04585 




















Withôt, ôx, and ôzinTable (1), Sandwas calculated and 
wave celerity C = a was calculated. The calculation results 


are presented in Table (2). 


Tabel.2: The value o and, and wave celerity C = = 





















































ôx ôZ o 

T ot St E 
(sec.) (m/sec) (m/sec) (m/sec) 
6 8,0677 8,05695 8,0677 
7 8,54589 8,5345 8,54589 

8 8,86229 8,85049 8,86229 
9 9,08074 9,06864 9,08074 
10 9,23739 9,22508 9,23739 
11 9,35337 9,34091 9,35337 
12 9,44158 9,429 9,44158 
13 9,51022 9,49754 9,51022 
14 9,56465 9,5519 9,56465 
15 9,60854 9,59574 9,60854 

It is interesting that = ~ = C.This correlation does not 


only occur for dispersion equations (11). If (11) is 
changed, it becomes: 

y?o? = gk tanh kh 
Whereyis a positive number greater than one, wavelength 
resulted will be shorter and the relation of z = z =C is 
obtained. 
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2.3. A function of two variables f (x, t) = cos kx cos ot 


The form of Taylor Serieswith two variables with 
variables(x, t), to ease the writing, it can be written: 


f(t + 6t,x + 6x) = ftx) +5, +52 +53 -t Sy 
..(12) 
S is the first differential term, szis the second differential 
term,and so on 
Next, the optimization equation is made: 


S2+5S3 


< ¢...(13) 








Functionf (x,t) = coskxcosot,is substituted tos,, S3, 
and sto (13) and made at a characteristic point 





wherecos kx = sin kx = cos øt = sin øt. The polynomial 
equation fordxis: 


Co + ôx + €,6x? + c,6x3 = 0....(14) 


_ ôt? ôt? i 
Co = 0 a a E 
ôt? 

c =-|o z totte k 
k2 
= (1 — oôt) — 
C2 = ( oO IZ 

k? 
tmr 


The equation can be written into an equation forôt, 
However, in this study, the equation is made with input 
ôtto calculatedx, wheredtis calculated with (8). Table (3) 
shows the calculation results of wave number kcalculated 
by the dispersion equation of the linear wave theory (11), 
with water depth ofh =10 m. 


Table.3 :The results for the calculation ofôtandôx with 



































(14) 

T ôt ôx 

(sec.) (sec.) (m) 
6 0,00191 | 0,04628 
7 0,00223 | 0,05719 
8 0,00255 | 0,06778 
9 0,00287 | 0,07813 
10 0,00319 | 0,08831 
11 0,0035 | 0,09836 
12 0,00382 | 0,10831 
13 0,00414 | 0,11819 
14 0,00446 | 0,12801 
15 0,00478 | 0,13779 
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Withdtandédxin Table (3), tig calculated with the 


calculation in Table (4). 


ôx 


; OX / 5¢ 
Table.4: The value o al 
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T = Car 

t k u ; 

(sec) (m/sec) (m/sec) C 

6 24,2138 | 8,0677 | 3,00133 
7 25,649 | 8,54589 | 3,00133 
8 26,5987 | 8,86229 | 3,00133 
9 27,2543 | 9,08074 | 3,00133 
10 27,1245 | 9,23739 | 3,00133 
11 28,0726 | 9,35337 | 3,00133 
12 28,3373 | 9,44158 | 3,00133 
13 28,5433 | 9,51022 | 3,00133 
14 28,7067 | 9,56465 | 3,00133 
15 28,8384 | 9,60854 | 3,00133 








5x3 a3 
MaA 
6 0x3 
ôx? arf 
— §z ——— 
ty ax2az 
z? af 
+6x —— 
*D axdz? 
Substitution, 
f(x,z,t) = cos kx cos ot cosh k(h + z) 
tos;, Szandszand optimization equation done at 


characteristic points and in conditionscoshk(h + z) = 
sinh k(h + z), equations fordz was obtained, wheredtand 
xas input, dtwas calculated using (8) while dxwas 


calculated using (14), 

Co + C1 0Z + C2627 + C3627 = 0....(15) 

With elements ofc, c1, cpandc3in Table (6) 

The condition coshifok (h + z) = sinhifoik (h + z) can be 
obtained in deep water. However, it does not mean that the 
obtained equation only applies to deep waters, it also 


applies to This is considering the 
conservation law of the wave number (Hutahaean (2020): 


shallow waters. 


OK(h+z) _ 


In contrast to the results of separatecalculations, using 
equations derived from equationsf(x,t)it was obtained 


that = = 3.00133 C,this fits the criteria of Courant (1928) 


ôx 
that — = 3C. 


1.4. A function with threevariables 
f(x, z, t) = cos kx cos ot coshk(h + z) 


The Taylor series for a function with three variables up to 
isf (t + ôt, x + 6x,z+ 6z) = 


the third 
f(t,x,z) + sı + 


derivative 


S2 + S3 


with s,, S2, and s3in Table (5) 


Table.5: Element S4, S,andsz 









































Sı S2 S3 
of St? a2 f St a'f 
ôt — a nearer 
ot 2 ot? 6 ot 
af af Se of 
+ôx — 2o Aa 
ax toet Sea. Foe aay 
af a2f ôt? af 
+6z— a sz—— 
Zaz TOLOA ae tz raz 
5x? a? bx? 93 
pa ig 
2 Ox? 2 ôtðx? 
a?f af 
+ôxôZ a +5t6x0z A 
ôz? ð? z? 8? 
B pa 
2 dz? 2 dtdz? 
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a = On ase (16) 
Table.6: Elementsof Co, C1, C2, and C3 
Co Cy C2 C3 
eodt —ek k? k? 
2 6 
+ekôx —okôt ok? 
——— ôt 
2 
j ôt? —k?8x k? 5 
o 5 5 x 
+okôtôx z, ôt? 
—g?k — 
2 
ôx? +ok?°ôtôx 
me ae 
2 
ot? ôx? 
+0? — —k? — 
gr 2 
ôt? 
+o7k —— ôx 
2 
x2 
+ok?ôt — 
ôx? 
+k? — 
6 
Table (7) shows the calculation resultsof ôt, 


ôx,andôzwherekwas calculated by the dispersion equation 
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of the linear wave theory (11), with the water depth 
ofh =10 m. 


Table.7: The calculation results of ôt, 6xandéz 



































T ôt ôx Oz 

6 0,00191 0,04628 0,13914 
7 0,00223 0,05719 0,17195 
8 0,00255 0,06778 0,20379 
9 0,00287 0,07813 0,23491 
10 0,00319 0,08831 0,26551 
11 0,0035 0,09836 0,29573 
12 0,00382 0,10831 0,32566 
13 0,00414 0,11819 0,35536 
14 0,00446 0,12801 0,38489 
15 0,00478 0,13779 0,41427 




















Withôt, ôxandôzin Table (7), and was calculated with 


the results presented in Table (8) 


Table.8: The calculation results o ~and=andC 















































ôx Oz o 

T St a | eok 
(sec) (m/sec) (m/sec) (m/sec) 
6 24,2138 72,8026 8,0677 
7 25,649 71,1177 8,54589 
8 26,5987 79,9729 8,86229 
9 27,2543 81,9442 9,08074 
10 27,7245 83,3578 9,23739 
11 28,0726 84,4044 9,35337 
12 28,3373 85,2004 9,44158 
13 28,5433 85,8197 9,51022 
14 28,7067 86,311 9,56465 
15 28,8384 86,7071 9,60854 








+1, OX 6z š ôx 6 ôz/ 6 bz 
With—and— dan Cin Table (8) —%, —and—was 
st. St c c ôx 


calculated with the results presented in Table (9). 


dx/ 6z/ 
Table.9: The calculation results a — Stand 











6z 
ôx 
T 
On), OF) x, bz 
(sec) C C ôx 
6 3,00133 | 9,02395 | 3,00665 
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7 3,00133 | 9,02395 | 3,00665 
8 3,00133 | 9,02395 | 3,00665 
9 3,00133 | 9,02395 | 3,00665 
10 3,00133 | 9,02395 | 3,00665 
11 3,00133 | 9,02395 | 3,00665 
12 3,00133 | 9,02395 | 3,00665 
13 3,00133 | 9,02395 | 3,00665 
14 3,00133 | 9,02395 | 3,00665 
15 3,00133 | 9,02395 | 3,00665 




















Referring to the calculation resultsin Table (9) relations 
can be formulated: 


ôx = 


Tk 


o 


y ôt.....(17) 
6z = a OF seins (18) 


With separate claculation as a function of single variable, 


8 8 A ‘ 
= = x =C was obtained, ory = 1 whereas with the 


k 7 6. 6. ; : 
simultaneous calculation = = 3C and = 9Cis obtained 
. 6 
ory = 3, all are related to wave celerity C.Thus,to make = 


closer to horizontal velocity u and ~ closer to vertical 
velocity w, Weighted Taylor series O(6+)on sinusoidal 
function f (x, t)should be in the form, 


f(x + 6x,t + 6t) = f(x,t) + yet + ôx Z a9) 


Where the total acceleration obtained, 


DF 4, OF y OF 
Aey tuž... (20) 


Meanwhile, for the functionf (x, z, t), the form of Taylor 
series 0(6")is, 
f(x + 6x,z+6z,t+6t) = 
OP gp OF». git 
f(,z,t) + y*6ot oe yox aT ôz zD 


With total acceleration, 


Df _ 2f of of 
a deer + yun + wzo 
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IW. WEIGHTED CONTINUITY EQUATION, 
EULER’S MOMENTUM EQUATION,AND 
KFSBC 


3.1. Weighted Continuity Equation. 


j w(x,z + dz) 
u(x, z) u(x + ôx, z) 
dbz —> —> 
lates 
¢+______> 
ôx 


Fig.1: Control Volume for the Continuity 


Equation Formulation 
At t = t, thus (22) can be written, 


a. For constantz 
f(x + 6x, z,t) = f(x,z,t) + yôx oF (23) 


b. For constantx 


f(z + 62z,t) = f(x,z,t) + 5z Z.....(24) 


The law of conservation of mass for the volume of a 
constant control volume (Fig. (1) and for incompressible 
flow, 

I—-0=0 

I = pu(x,z,t) 6z + pw(x,z,t) 6x 


0 = ( t) + oxo) 
= p\|u(x,z,t) yox a Z 


+ ( (x,z,t) + 6 Wa 
p|w(x,z, Zz)? 


Subtraction and equation are divided bypôxôz, weighted 
continuity eqation is obtained, 
ðu 


ðw 
Yota } ee (25) 


3.2.Weighted Euler’s Momentum Equation 


Using (22), weighted Euler’s Momentum Equation in the 
direction of the direction of horizontal axis-xand in the 
vertical axis-zare 


y u rTu pi signe scree (26) 
2 dw dw Want 
ra ore a ead | (27) 


3.3. Weighted KFSBC. 
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The known KFSBC (Dean (1991) is, 

Wy = ou + Uy - 

w,is the water particle velocity on the surface which is the 
total velocity of the water level elevationy(x,t) = 
coskxcosot, while theweighted total acceleration of 
water level elevation with(20) is 

Dn on ôn 

dt Vat t ax 

Thus, wigthed KFSBC is, 





ð ð 
Wy =Y SE + Uy se. --(28) 


THE APPLICATION INTIME SERIES WATER 
WAVE MODELING 


In this section, the governing equations for time series 
modeling water wavesare formulated. The governing 
equations consists of two equations, they are the water 
surface elevation equation and the particle velocity 
equation. The variable of particle velocity in this equation 
is the depth-averaged velocity. 

a. Water surface elevation equation 

The Continuity equation (25)is multiplied by dzand 
integrated with water depth. . Integration of the first term 
is completed with the Leibniz integral (Protter, Murray, 
Morrey, Charles, 1985).KFSBC and bottom boundary 
condition were calculated, 


a a a 
y= fi u dz- (y — lu, +y? = 0...29) 


The integration of the left-hand first term is solved by 
using the particle velocity equation for the solution of 
Laplace’s equation (2). From (2), the relation of the 
direction of horizontal axis of particle velocity at an 
elevation z to the horizontal velocity at elevation y is 


_ coshk (h +z) 
u = cosh k(h+n) “a 
Left hand integration (29) becomes, 


f F -f coshk (h + z) 
Bi on _, coshk(h + n) Z Un 


Integration is completed using 7 = “and defined by H = 


h +î and calculated in deep water depth 
wheretanh k (h +n) = 1, 
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| udx = 

i k 


Conservationlaw of the wave number (Hutahaean (2020) 
is, 


dk (h + $) 
dx ~ 
or 
A Ay 
k(n +5) = ko (ho +2) 


In deep water tanh Ky (ho + 2) = 1whereky (ho + £) = 


Or, a relation is obtained 


On On 





Tog 


The final result of integration is, 


| d U,H 
udx = Dr 


-h 
Substitute to (29), 


on 1 OuyH 


a Ge ox T 





@-1), an 
7 Un zye 680) 


As mentioned earlier, the modeling uses depth-averaged 
velocity. The horizontal depth average velocity Uis 


defined as the particle velocity at the leevationz = 
Zobelow the SWL, wherezois a negative number. From (2): 


Un cosh k H 
U coshk(h +z) 





Is defined: 


cosh kH 
= Cosh k(h+z9) 
Thus, the relation of horizontal surface velocity with 
horizontaldepth-averaged velocityis: 


Up = AU (32) 


Substitute to (30) 
on _ a ðUH (yv-1) 





an 
ap a ae 7 aU —......(33) 


Zois calculated by the following equation, 


1 Als 


— dz=1 
mla Z 
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From (2) and the definition ofdepth-averaged velocity, 


_ coshk (h +z) 
s cosh k(h + Zo) 


The characteristic ofzois, 
1/2 coshk (h +z) 

— Í ——— ~ dzU=1 

UH J_, coshk(h + zo) 

From this equation the equation for Zois formulated: 


kH coshk(h + zo) — sinh kH = 0...(34) 


The calculation ofain (31)and in the calculation ofzo (34), 
in the deep water depth,requires deep water depth depth 
value hg. wave number kọis calculated bydeep-water 
weighted linear wave dispersion equation, 





292 
Ko = id 
g 
As deepwater depth: 
on Ao 
Oe ko 2 


Deep water depth is used to calculatea.Considering 
conservation law of the wave number, the value of qis 
constant. 


b. Horizontal velocity equation 


Weighted horizontal surface momentum equation 
(Hutahaean (2021), is, 
du, 1d dn 
n — 
2 at + Dax (yunun + WyWr) = —g ax 
By substituting surface velocity with (32) andby 
neglecting convective acceleration, 
I = --451_35) 


dat — ay ax" 


c. Model Results 


The Finite Difference Method for spatial differentials uses 
numerical solutions, while time differentials are solved by 
the predictor-corrector method for numerical integration 
(Hutahaean, 2019). The time step dtwas determined with 
(8), using €= 0.005, while the grid size ofdxwas 
calculated with (17).The model executiom was made using 
weighting coefficient y = 3.0, deep water coefficient 0 = 
2.0, tanh 6r = 0.99999. 


As the first case, the model was done on a channel with a 
constant depth h = hy = 14 m. In the channel there are 
sinusoidal waves with wave periodT = 8 sec. and wave 
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amplitude A = 1.0 m. The model results are presented in 
Fig. (2). Fig (2) shows that the model can simulate well the 
short waves with large amplitudes. 


1.5 
1 
0.5 
0 
-0.5 


n (m) & U (m/sec) 


-1 


-1.5 
0 50 100 150 200 


x (m) 





surface elevation vel. U 





Fig.2: Model Results on Flat Bottom 


In the next case, the model was made on a sloping bottom 
13 

E 

is họ = 14 m, while upstream water depth is 1.0 m. The 


incoming wave of sinusoidal wave with the wave period 
ofT = 8 sec. and wave amplitude ofA = 0.8 m. The model 
results are presented in Fig.(3). 


; dh 
with a bottom slope Ti Downstream water depth 


3 
25 
2 
1.5 
1 
0:5 
0 
-0.5 
-1 
alk5: 


n (m) & U (m/sec) 


0 50 100 150 200 


x (m) 








surface elevation vel. (U) 


Fig.3: Model Results on Sloping Bottom. 


The model results show that initially shoaling occurred, 
then the waves became unstable at waterdepth of 4 m and 
then the breaking peak occurs at a water depth of 2.60 m. 


V. CONCLUSION 


Some conclusions are drawn from this study. The first is 
that the application of the Taylor series to the sinusoidal 
wave equation should use time step and spatial intervals 
correlated with phase speed. Thus, it can be correlated 
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with the water particle velocity, a weighting coefficient 
must be obtained. The Taylor series with the weighting 
coefficient is hereinafter referred to as the weighted Taylor 
series that only uses the first derivative. 


The formulation of hydrodynamic equations with the 
weighted Taylor series produces equations with the 
weighting coefficient, including the weighted continuity 
equation, the weighted Euler’s momentum equation,and 
the weighted KFSBC. 


The next conclusion is that by using the weighted 
equations, the time series of the wave equation is obtained 
to simulate a shortwave where short wavelengths are 
produced and there is a breaking phenomenon. 


The determination of the time step and gridsize in 
numerical modeling using the Finite difference method can 
use the equations formulated in this study. 
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